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PREFACE 


This report details the modeling phase of the project "Modeling and 
Control of Transonic Cryogenic Tunnel" sponsored by NASA/Langley Research 
Center under grant NSG 1503 which was performed during the period ending 
September 1979, The scope of this report is confined to model synthesis 
and its validation. Preliminary findings on this modeling activity w^re 
reported at the First International Svinposium on Cryogenic Wind Tunnels 
which was held at University of Southampton, United Kingdom, in .April 1979 
(ref. 2) . 
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SYNTHESIS OF A CONTROL MODEL FOR A LIQUID 
NITROGEN COOLED, CLOSED CIRCUIT, CRYOGENIC, NITROGEN 
WIND TUNNEL AND ITS VALIDATION 


By 

S. Balakrishna^ 

INTRODUCTION 

Demand for full-scale Reynolds number testing in modern transonic wind- 
tunnel test programs has resulted in considerable review of various methods 
of realising such high Reynolds numbers without other serious penalties. The 
most promising amongst many techniques has been the concept of operating the 
test gas of a tunnel at cryogenic temperatures (ref. 1) . Due to advances in 
the cryogenic engineering of economical production of cryogens, materials 
compatible for cryogenic operation, and instrumentation, worldwide efforts at 
development of high Reynolds facilities using closed circuit cryogenic pres- 
sure tunnel development can be noted (ref. 2) . All this activity was spear- 
headed by the proof of concept operation of the 0.3-m transonic cryogenic 
test facility at NASA/Langley Research Center (LaRC] (ref. 3). Operational 
experience with this facility has shown that cryogenic operation of its tunnel 
gases is a practical method for obtaining high Reynolds number flows without 
any dynamic pressure penalty and with savings in fan power (ref. 1) . 

Confidence in a wind tunnel is related to the flow quality and to the 
stability of the cunnel variables during the test period. Initial tunnel - 
running experience on the 0.3-m transonic cryogenic test facility has indicated 
that, though the tunnel process was basically stable, continuous drifting of 
the tunnel variables occurs because of mass and energy disturbances caused by 
unmatched control inputs. The instability is further compounded by a consider- 
able degree of cross coupling between the tunnel variables and the controls. 
Manual regulation of the tunnel variables, though feasible, proves to be rather 
inefficient in terms of liquid nitrogen consumption. Further, transition of 


^ Research Assistant Professor, Old Dominion University Research Foundation, 
P.O. Box 6369, Norfolk, Virginia 23508. (On sabbatical leave from National 
Aeronautical Laboratory, Bangalore - 560 017, India.) 


the tunnel conditions from one state to another has also proven to be uneconom- 
ical, due to long periods necessary to realize acceptable stability of tunnel 
conditions. Obviously closed-loop regulation and control of the cryogenic 
tunnel are the only solutions to the control problem. 

For these reasons, studies have been undertaken to develop a control com- 
patible mathematical model of a typical closed circuit cryogenic tunnel. No 
literature on such a model exists since the concept of cryogenic operation of 
tunnels is a recent one. In this document, the details of the efforts to 
synthesize a control -compatible multivariable model of a liquid nitrogen cooled, 
gaseous nitrogen operated, closed circuit, cryogenic pressure tunnel are 
presented. The synthesized model has been transformed into a real-time cryogenic 
tunnel simulator, and this model is validated by comparing the model responses 
to the actual tunnel responses of the 0.3-m transonic cryogenic tunnel, using 
the quasi-steady-state and the transient responses of the model and the tunnel. 
The global nature of the simple, explicit, lumped multivariable model of a 
closed circuit cryogenic tunnel is demonstrated. 
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NONENCUTURE 


A area of cross section Ctn^) 

C velocity of sound (m/ ) 

Cy specific heat at constant volume (J/kg K) 

Cp specific heat at constant pressure (J/kg K) 

specific heat of tunnel metal (J/kg K) 

Cj valve coefficient 

d derivative operator 

E energy (J) 

g gravitational constant 

H,h enthalpy (J/kg) 

i,I integer iT'^.ices 

J mechanical equivalent of heat 

k,K constants 

i length (m) 

LN'2 liquid nitrogen 

GN2 gaseous nitrogen 

• 

m mass flow (kg/sec) 

M Mach number 

N fan speed 

P tunnel pressure (atm) 

Q heat (J) 

q dynamic pressure (kg/m^) 

r fan pressure ratio 

R universal gas constant 

Re Reynolds number 


j tharmal resistance [(J/m^K sec)'M 

5 Laplace operator 

T gas temperature (K) 

t time, time constant (sec) 

U internal energy (J) 

u specific internal energy or velocity or m/sec) 

V volume (m^) 

W mass (kg) 

X distance (m) 

y,Y heat transfer coefficient (J/m^K sec) 

Z thermal impedance 

a cooling capacity of gas bleed (J/kg) 

6 cooling capacity of liquid nitrogen (J/kg) 

Y ratio of specific heats 

p density of gas (kg/m^) 

A increment 

J thermal mass (kg) 

u viscosity 

T transport delay (sec) 

n efficiency 

• dot over a quantity refers to ^ of that quantity 

Subscripts 
F,f fan 
g gas 

k kinetic 

L liquid nitrogen 
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Subscripts (concluded) 
m metal 

t tunnel 

s static 
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PRINCIPLE OF OPERATION OF SINGLE- RETURN, 
CLOSED CIRCUIT WIND TUNNELS 


The wind tunnel is a device intended for studying the behavior of scaled 
models in fluid flow. A single-return, closed circuit wind tunnel is an 
efficient and aerodynamical ly designed endless duct made up of four straight 
segments which are joined by four corners. The wind-tunnel resident fluid, 
usually air, is moved axially by means of a fan to achieve desired flow veloc- 
ities at a segment of the tunnel called the test section, wherein the test 
models are located. The geometrical design of a typical single-return, closed 
circuit tunnel is carried out such that the desired range of flow velocities 
is achieved at the test section with minimal fan power requirements and the 
flow possesses the desired flow qualities (ref. 4) . Such tunnels typically 
consist of a settling chamber, a screen section with a heat exchanger to reduce 
turbulence and remove heat, a contraction which allows smooth and efficient 
acceleration of the test gas to desired flow velocities, a test section with 
a plenum to obtain adequate control of wall interference, a return leg diffuser, 
the first and the second leg comers with turning vanes, a fan and fan nacelle 
•section with fan motor located either inside or outside the tunnel, a forward 
leg diffuser, and the third and fourth leg corners with turning vanes leading 
into the settling chamber. The closed circuit tunnels are either open to 
atmosphere or are totally enclosed, and these are called atmospheric or pressure 
tunnels respectively. 

In order to accurately extrapolate the wind-tunnel generated, scale-model 
data to their full-scale counterparts, it is necessary to maintain flow simi- 
larity conditions in the wind-tunnel flow (ref. 4) . The important flow simi- 
larity conditions for stationary test models are Reynolds number, the ratio of 
inertial and viscous forces, and Mach number, the ratio of inertial and elastic 
forces. Dynamic pressure is another flow parameter of model strength and 
structural significance. These flow parameters Re, M, q, are coupled non- 
linear functions of the measurable test section flow variables, viz static or 
stagnation pressure, static or stagnation temperature, and the mass flow rate. 
Under isentropic flow conditions, the tunnel test section flow parameters and 
the flow variables are related as 
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The three flow parameters Re, M, and q are functions of three tunnel flow 
variables P^, and m. To maintain flow similarity condition^' in model 

study, a wind tunnel must allow independent control of the three fi .-; parameters, 
which in turn can be achieved only if the three tunnel variables P^, T^ and 
III can be controlled independently. 

Historically, in wind-tunnel testing, flow similarity has rarely been 
achieved in respect to Reynolds number. A study of the flow parameter control- 
lability in various types of tunnels can now be reviewed. .Ambient temperature 
atmospheric tunnels have an uncontrollable P and T which are always near 
atmospheric. In such tunnels only m and hence Mach number can be independ- 
ently controlled. Naturally, such tunnels are associated with uncontrollable 
Re and q. In the case of ambient temperature, closed circuit pressure 
tunnels, sometimes called the variable density tunnels, the test mediitn temper- 
ature is near atmospheric. However the tunnel pressure P can be controlled 
by adjusting the mass of fluid in the system, which is in addition to mass flow 
and Mach number control. Thus the variable density tunnels allow independent 
control of m and P and hence M and q, but have a dependent Reynolds 
number. A high Re can be had, but is associated with considerable increase 
in flow dynamic pressure q and its accompanying structural problems. It is 


in this context that the cryogenic operation of the test gas of a closed circuit 
pressure tunnel has become relevant. Such a tunnel allows independent control 
of all the variables m, P and T and hence independent control of M, q 
and Re . 

The test medium gas temperature control problem is general to all closed 
circuit, fan-driven tunnels, since operation of the fan to create desired 
flow is associated with fan compression and wall friction heating which tend to 
increase the gas temperature continuously. In the case of ambient temperature 
tunnels, this heat is removed by either a heat exchanger of the water tube t\'pe 
or by mixing atmospheric air. In either case the temperature control is crude 
and at best can provide cooling down to ambient. Further, the heat removal 
schemes are bulky and do not provide an adequate range of temperature control to 
influence Reynolds number. 
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CLOSED CIRCUIT, LIQUID NITROGEN COOLED, GASEOUS 
NITROGEN OPERATED, CRYOGENIC PRESSURE TUNNELS 

The cryogenic wind-tunnel concept (ref. 1), though proposed many times, 
has become an engineering reality because of advances in cryogenic engineering 
and material. The concept allows efficient cooling down or regulation of the 
tunnel gas temperature at any fan power and involves injecting a liquified form 
of the test gas into the tunnel. Usually, gaseous nitrogen (GN 2 ) is the test 
medium and liquid nitrogen (LN 2 ) is the coolant. This cooling mechanism is 
fast and efficient, and cooling occurs down to near vapor phase temperature. 

In case of nitrogen, this temperature range is from about '8 K to 350 K at 1 
atm. In contrast, the ambient temperature tunnels operate at a few degrees around 
the ambient 325 K. The wide range of temperature control in the ciyogenic tunnel 
allows a large increase in Reynolds number without any dynamic pressure penalty. 
Since the viscosity of nitrogen gas varies as T®*^, the Reynolds number 
increases as a function of T"^***. This is equivalent to a si.xfold increase 
in Reynolds number when the gas is cooled from 500 K to 100 K. 

Under steady-state operating conditions, the cooling induced by liquid 
nitrogen injection equals the fan power, thereby maintaining the temperature 
and keeping the tunnel pressure constant. Addition of liquid nitrogen into 
a tunnel not only results in cooling, but is also associated with tunnel gas 
mass increase. In order to maintain the tunnel resident gas mass, it is 
necessaxy to remove an equal amount of warm tunnel gases, which is associated 
with removal of enthalpy. This complex interaction of fan-induced heating, 
liquid nitrogen contributed mass increase and enthalpy decrease, and gaseous 
nitrogen bleed-induced mass decrease and enthalpy decrease is the cryogenic 
tunnel control problem. This can be analysed by studying the thex'modynamics 
of the cryogenic tunnel and its control inputs, which are intimately related 
to the thermophysical properties of liquid and gaseous nitrogen. 

Thermophysical Propex'ties of Nitrogen 

The physical and thermal properties of element nitrogen have been a subject 
of e.xhaustive studies over wide ranges of temperature, pressure, and both 
liquid and gaseous phase. These properties have been well documented (I'ef. 5). 


These studies indicate that under static conditions, nitrogen deviates con- 
siderably from perfect gas behavior in many respects. However, under 
isentropic expansion or contraction, nitrogen gas has been proven to behave 
near perfectly (ref. 6), both thermally and calorifically. 

Real gas static properties, particularly density p and specific heat at 
constant volume C^, are of special importance to control analysis. These 
real gas variables p and are nonlinear functions of temperature and 

pressure, and the nonlinearities are dominant near the liquid-gas phase 
boundary. A simple function fit has been made for p and using data 
from reference 5. 


P = 12.18 Y 
or 

1 + 250 

T2J 

mol/liter Sat £ T £ 400 K | 

p = 338.9 ^ 

1 + 250 -^1 

T^J 

kg/m^ 1 £ P £ 8 atm I 


Similarly 


C » 20.80 fl + 2S0 J/mol K Sat < T < 400 K 

V \ jZ J 


1 < P < 8 atm 


(3) 


(4) 


Figures 1 and 2 illustrate real gas properties p and respectively, as 

functions of T and P. Nitrogen in its liquid form is associated with an 
enthalpy hj^ which is a function of storage pressure and temperature. When a 
unit mass of liquid nitrogen stored at I atm pressure and saturation temperature 
is taken to a higher pressure P at temperature T, the heat absorbed by the 
unit mass of liquid nitrogen is the cooling capacity of liquid nitrogen. This 
enthalpy change has been estimated by Kilgore et al. (ref. 7) using the temperature 
entropy diagram. The route of this enthalpy change is a first isentropic 
compression of liquid nitrogen to the final pressure P, and subsequent warming 
and evaporation from colder saturation to higher T isobarically. From ref- 
erence 7, we have 


s « 


cooling capacity 
- - [bqT + Bi + ^] kJ/kg 

where 

Bo « 1.0379 - 3.91? x 10"3 p 

Bx » 121.25 + 2.157 P 

B2 = 66.85 - 391.2 P 

P = Pressure in atm 

T * Temperature in K 


This can be appro.x'imated to 


S “ h^^ - CpT kJ/kg 

Cooling capacity S is a negative number, where 
h, * Latent heat of vaporization - C T 

L P Sat 

C = Specific heat at constant pressure = 1.04 kJ/kg 
P 


( 5 ) 


A similar analysis of the cooling capacity of gaseous nitrogen taken from pres- 
sure P and temperature T to 1 atm and higher ambient temperature can be 
estimated on the temperature- entropy diagram. This can be approximated to 


a == CCp - C^) T kJ/kg (6) 

where 

a = cooling capacity of gaseous nitrogen 
Cp, Cy = specific heats of nitrogen 
T = Temperature, K 

the cooling capacity ct is a positive number. Cooling capacities, 3 of 
liquid nitrogen and a of gas bleed are shoivii in figures 3 and 4, respectively, 
as functions of T and P. 
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Energy State Diagram 


The internal energy associated with a finite volume of gas at a uniform 
pressure P, and temperature T is all energy associated with electrons 
at the molecular level and is a function of temperature as predicted by the 
third law of thermodynamics. At absolute zero, the internal energy and gas 
specific heat tend to zero. This can be e.xpressed using the classical 
equation 

U C dT 

0 

where 

C„ * ^ at constant volume 
V at 


Density of nitrogen gas, as discussed previously, can be used to estimate mass 
of gas W , for a given volume of gas V at pressure P and temperature T 

•S 


The internal energy associated with this volume of gas is 


W = 338 . 90 ^ I 1 + 250 

§ 


U = W C dT 

si V 

T 

Assuming r C^dT = C^T we have 

*o 

U = 255.49 P^l + 250 ^1 + 250 kJ/m^ 


r- 1 


( 8 ) 


A plot of gas internal energy U vs. its mass is shoivn in figure 3, and 

it illustrates the mass-energy relationship for a unit volume of 1 m^. The 
loci of constant temperature on the diagram are the lines of slope C^T radi- 
ating from the origin, and these are the isothermal lines. The loci of 
constant pressure are lines which generally run parallel to x-axis and are. the 
isobars. If the gas were truly perfect statically, the isobars would be truly 
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p&ralldl to the mass axis. The real gas isobars liave positive slope. The 
liquid-gas phase boundary runs near the 100 K isothermal line, and its true 
position has been ignored. 

Consider a representative point A on the diagram, corresponding to a 
finite mass of gas at a know pressure and temperature in a 1-m^ volume. If 
any extra gas mass or energy is added to this system, the point A appro- 
priately shifts graphically to account for the new mass and energy state. In 
the case of the cryogenic tunnel, the three control inputs and the thermal 
leakage, either from the ambient or the tunnel metal shell change the energy 
state of the system by moving it to new pressures and temperatures. Addition 
of liquid nitrogen constitutes a reduction of energy by negative enthalpy 
mj^hj^ and increase of mass by mj^. Bleeding of gas from the system similarly 
represents a reduction of energy by negative enthalpy m h , a mass decrease 
by m . Heat release from fan is a true energy input Qp and wall released 

heat -Q^. All these inputs can be represented as mass-energy vectors, and 
the sum of these control vectors provides the final state of system. This 
analysis can be carried out graphically by vector manipulation as shorn in 
figure 5. 

The same analysis can be carried out analytically, by studying the open 
thermodynamic system. 


Thermodynamic Model of a Closed 
Circuit, Cryogenic Pressure Tunnel 

The cryogenic pressure tunnel can be considered as a thermally autonomous 
pressure vessel of volume V containing nitrogen gas at a uniform pressure P 
and temperature T. The vessel is assumed to be ideally insulated, insulation 
being external to the metal wall. This vessel is assumed to be opened to a 
liquid nitrogen source and a gas bleed valve which allow transfer of mass and 
associated enthalpy. The vessel is also assumed to be opened to a fan shaft 
which creates necessary fluid flow and delivers the fan energy. Such a system 
is shown in schematic in figui'e o. The liquid nitrogen flow rate is taken to 
be -and occurs at a liquid pressure of P^. The gas valve allows flow 

of m out of the tunnel. 

S 
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Assuming uniform mixing and ignoring work done, potential energy and 
kinetic energy and invoking first law of thermodynamics we have 




then 


where 


h„ » + P V 

S g g g 


* "g^ * 


-Q^ = rate of heat release from tunnel wall to gas, J/sec 

Qp = rate of heat release from fan, J/sec 

U = internal energy = W^u, J 
u = specific internal energy 

j dW , 

“Q,. + Q„ + m, h, ® m h + -jr- (W u) = m h + u + W 
n L L g g dt V g / g g dt g dt 

• • *. /• *\. du dT 

"Q* + Qp + m.h, - m h = u(m, - m ) + W, — • ;r- 

^t L L g g \ L g' g dT dt 


which can be simplified to 

"l ("l - - "g (i>g - - Q, * Qp = Yv § 

For a perfect gas we have u = C T and h * C T; hence 

^ V g p 


dT 


"l (^L - V) - "g (S - <=v) * «F - Vv H 


( 9 ) 


The mass and energy crossing the boundary of the tunnel react with tunnel resident 
gas and change its state. This interaction, as detailed previously, of the 
control vectors can be shown on the energy state diagram as the charge, bleed 
and fan mass-energy rate vectors. In order to run the tunnel at a steady gas 
state, the Mach dependent Qp needs to be cancelled by sum of negative enthalpies 
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from charge and bleed, which in themselves should create no mass change. By 
manipulating the magnitudes of these control vectors, either a balance condi- 
tion can be obtained or the tunnel state moved along any desired path. 

The details of vector manipulation are shown in figures ^ to 11. Figure 
7 shows conditions needed to maintain the steady tunnel state. The control 
input relations are 


'"l • "g ' - t » 0 . mL(hj^ - CpT) - 

In figure 8, showing the technique of increasing gas temperature at constant 
pressure, the control input relationships are 


m 


L 


and 


T 

g P 


JL T 
W 



P s 0 


In figure 9, the technique of gas temperature reduction at constant pressure 
is demonstrated. 


Aj = p . 0 


\'’l - Qp ■ ° 


T = 



In figures 8 and 9, has been taken to be zero. In figure 10, the manip- 
ulation of control vectors to realise a pressure increase at constant temper- 
ature is illustrated. The control input isolations are 


T - » 0 . mj^(hL - Cj) + Qp » 0 , P = m^C^T 

Similarly figure 11 illustrates control vector manipulation to achieve a pressure 
reduction at constant temperature. 


t = \ = 0 



P » - mgC/ 


IS 


Role of Metal-Stored Energy in Cryogenic Tunnels 


Cryogenic tunnel engineering design demands insulation of the cold test 
gas from the ambient to prevent heat gains and the associated energy waste. 
Since these tunnels are required to operate over a considerable pressure 
range, typically one to eight atm, the insulator shell must possess the ability 
to withstand internal pressure without leakage of gases. Most insulators are 
poor structural material, and storage of high-pressure compressible fluids in 
the shell demands perfect safety and conformity to appropriate pressure vessel 
code. These requirements dictate use of cryogenic structural metals for the 
pressure vessel integrity. However the location of the insulation on this 
vessel can be either internal or e.xternal. When the insulation used is 
external to the tunnel shell, the structural metal is e.xposed to cold gases 
and thermally interacts with the gas by releasing heat. Alternately, the 
insulation can be internal to the tunnel, thereby isolating the metal from 
cold gases. In the latter case, the need for very good surface finish on the 
insulator to obviate problems of boundary layer growth demands use of thin, 
highly finished, metal panels. This is particularly true of the test section 
and contraction wherein good aerodynamic contouring is a must. Thus, engi- 
neering design of cryogenic tunnels dictates use of a considerable amount cf 
metal within the insulation shell. 

When the insulator behaves ideally, the cryogenic tunnel can be looked 
upon as an autonomous thermodynamic system. The metal -stored energy then 
becomes a part of the tunnel thermal inertia and tunnel internal energy. 

The tunnel energy is now in the form of gas molecule motion and lattice 
vibration of metal molecules. Any gross temperature difference between gas 
and metal results in heat transfer between the two at the gas -wall boundary, 
and under steady-state conditions the metal and gas are at the same tempera- 
ture. The magnitude of metal-stored energy can be estimated by studying the 
enthalpy and specific heat of metals (ref. 8) 

Specific Heat of Metals 

The heat or energy stored in a molecule of solid metal at constant 
pressure is a function of temperature and can be defined as 
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where C « ^ J/raol, and C „ is the specific heat of metal at constant 

pressure P. The third law of thermodyTiamics predicts that for all phenomina 

which are in thermal equilibrium, the associated specific heat becomes zero at 

absolute zero temperature. For most engineering purposes, in the range of 

temperatures relevent to crvogenic tunnels, C * C = C . Specific heats 

mp mv m 

for ccnst int pressure and volume are almost equal. 

It is well known that metals with a face-centered lattice crystal structure 
are the best-suited metals for cryogenic operation. Aluminum alloy bObl-Tb 
is a well-proven cryogenic aluminum alloy. In the rest of the analysis, use 
of this alloy has been assumed. Enthalpy of aluminum as a function of temper- 
ature has been well researched and reviewed by Corrucini (ref. 9). The specific 
heat data is also available in the same reference for the tempex*ature range of 
cryogenic tunnel interest. Simple function fits have been made: 

C = 5.5 T - O.OOS T- J/kg K T < 350 K (10) 

m ^ ' 


H Cj^dT = 2.75 T- - 0.00266 T^ J/kg T < 350 K (11) 

0 

These are illustrated in figures 12 and 13. 


Cryogenic Tunnel Energy State Diagram 

Consider a closed circuit, cryogenic pressure tunnel made of aluminum 
which is externally insulated. The total energy contained within this system 
is the sum of gas internal energy and metal enthalpy: 


E = U . 


W C T + 

g V 



C dT 
m 


( 12 ) 


The energy state diagrams for two different, closed circuit, cryogenic, pres- 
sure tunnels, viz, the 0.3-ra transonic cryogenic tunnel (TCT) and the proposed 
National Transonic Facility (NTF) are now analyzed. The 0.3-m TCT is an 
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externally insulated tunnel of total volume 14.1 and 3200 kg of aluminum 

alloy 6061-T6 within the insulation shell. The NTF has an estimated volume 

of 5600 m^, an internally insulated, stainless steel pressure vessel, and an 

estimated 280,000 kg of aluminum liners turning vanes, contraction cone, 

test section and other components. Two sets of the energy state diagram, 

for both the 0.3-m TCT and the NTF are presented in figures 14 to 17. Figures 

14 and 15 correspond to U vs. and E vs. for the 0.3-m tunnel.. 

Figures 16 and 17 are U vs. W and E vs. W for the NTF. In the 

S S 

case of the 0.3-m tunnel, the maximum U is only 3 percent of E and, in 
the case of NTF, U is about 12 percent of E. In either case, the 
dominant energy content of metal is evident; hence the metal thermal inertia 
is likely to be the dominant factor which affects the pressure and tempera- 
ture dv-namics of the tunnel. The figures 14 to 17 also illustrate the 
control vector directions. It can be noted that a simple charge input in 
figures 14 and 16 shifts the tunnel states such that both temperature and 
pressure drop. However, in figures 15 and 17, the same charge input generates 
a negative temperature shift and a positive pressure shift. Thus the energy 

state diagram U vs. W provides the instantaneous behavior of tunnel 

S 

states and diagram E vs. W provides the final steady-state response. 

o 

As will be demonstrated later, this shift of response from U initially 
to E provides the characteristic signature of cryogenic tunnels. 


Metal to Gas Heat Transfer 

The mechanism of heat transfer in a tunnel is a surface phenomenon and is 
a complex function of the metal and gas properties such as thermal conductivity, 
specific heat, gas density, and gas viir.csitv. Assuming the metal to gas heat 
transfer to be a forced convection heat transfer and ignoring other modes, the 
following dimensional analysis can be made: 


Density of gas 

P 

ml’ 3 

Specific heat 

C 

P 

L^t’^T’l 

Thermal conductivity 

1/R 

MLt’^T’l 

Viscosity 

u 

.ML’U’l 

Velocity 

u 

Lt’l 
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Distance i 

Heat transfer coefficient y 


L 

Mt'^T’l 


T!ve heat transfer coefficient y, is now assumed to be a function of the above 
listed variabl'^s. Matching the dimensions, we have 


MT'^t"^ » A 0 Cp u u'" 

using arbitarary indicies £, b, £, £, For a dimensional match to be 

true 


a b 


^ /Ai 


(pu)^ (uc„r)1 


DR 

For the following relations of the gas density, viscosity, and thermal conductivity 

D « P_ ; u « m/f ; UocxO.9 . r® 7 * 0.98 
T 

and Prandtl number uCpR is assumed to be constant. Then 

pO.8 7~0.204 [t|0 . 8 
y 

dO-2 

assuming b to be 0.8 (ref. 10). A numerical estimation of the heat transfer 
coefficient is now made based on constants validated by Bartz in reference 10. 

Using equation (6) of Bartz 's work and using his notation. 


h = 


0-2 c \ / P V-e-i /a*V-9 


DO.2 \ prO.a /\c-^ 


a = 


T V 

- 

0.632 

0 .167 

1 w / 

2 T V 

0 

1 + 0.2M2'^ + j 


L + 0.2M2^) 


using constants of w = 0,6 and C = 0.026. This heat transfer coefficient 
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has been validated at the nozzle of a rocket. A number of values of were 
estimated using Pr « 0.73 and C* ■ x 61.5 for various pressures, 

temperatures and area ratios. This set of data was function fitted to obtain 
a simple heat transfer equation 

Yg » 71.27 PO-® t' 0-22 M(1 - 0.67M J/m^ K sec (13) 

This heat transfer coefficient Y^ provides the local forced convection heat 
transfer on a unit area of the tunnel metal-gas surface and is a function of 
local pressure, temperature, and Mach number. It can be used to establish 
the overall heat flow rate. 

The heat flow from wall to gas is not only a function of surface convec- 
tion based Y , but also a function of wall thermal conductivity, wall 
thermal capacity and wall leakage, if any. Consider a unit area of the tunnel 

i 

wall thickness, starting from the gas boundary to the external tunnel face, as 

shown in figure 18. This cross section consists of the tunnel metal wall and ; 

the insulation. The heat flow occurring under transient conditions can be 

obtained from the simple electrical analog. Here 1/R coiresponds to the 

heat transfer coefficient Yg associated with the turbulent boundary layer 

of the tunnel flow discussed previously. Further, corresponds to the j 

thermal conductivity of the metal wall thickness and corresponds to j 

thermal capacity of that layer. In view of the serial disposition of metal, j 

R is divided into a number of layers each of resistance R, and C„ is i 

divided into the same number of layers, each of capacity C. Similarly | 

1/R^ corresponds to thermal conductivity of the insulation, which tends to 

1 i 

zero, and Cj is the enthalpy in the insulation. The heat-flow dynamics ! 

can be now analyzed by evaluating the thermal admittance of the ladder network. j 

1 » ^ 

Z = Impedence = R + / 

CS + 1 

R + 1 

CS + 1 . 

R + 1 

CS + 1 

R + ... 
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The above ladder network has been simplified for five serial elements of 
metal ignoring insulator: 


iM 


3CS 


f f :3(Rg^s^ ^ ssiro-s- t- isrcs i 


siRci^s^ ^ :uRci“S‘ ^ :o(Rcis ^ ^ 


ii4) 


The time constant of one wall is _'5RC « Using this a pole ;oro analysis 

of the impedence function was made. The poles are 

dS 


Poles 


i-O.S * j O.ol 

- 59 

TT~ 

m m 

- 100 

R C 
m m 


R C 
m m 


i « v-l 


Zeros 


R 0 
m m 

- r.5 

R C 
m m 

- 45 


V'm 

- 90 


R C 
m m 


Ignoring the larger values of poles and :ero, a simpl ification is made to esti- 
mate Z. 

R C Ss 




(151 


Combining the wall and gas heat transfer equations, the heat flow in a unit 
area of tunnel wall on the basis of one-dimonsional he.it flow, and ideal 
insulation: 


:i 


( 16 ) 



Equation (16) provides a simple, lumped model of heat flow in each unit area of 
the tunnel surface. The total heat flow can be estimated by integrating the 
individual heat flows. This is fairly comple.x in view of the fact that Y » 1/R_, 

o o 

the dominant element of equation (16) , varies throughout the tunnel surface 

R C 

as detailed by equation (13) . It may be noted that lead time constant m m 

20 

is a small one compared to the denominator, and hence the heat flow is basically 
a first order dynamic. 

The assumption of one-dimensional heat flow only across the wall, and not 
along the surface, can be justified on the basis that the 0.3-m TCT is flanged, 
and these flanges occur wherever area is changing rapidly. The flanges carry 
insulation in the form of a pressure seal gasket and hence prevent transverse 
heat flow. In any such flanged segment the tunnel states P, T and M are 
uniform to a first ox'der appi'oxiraation . All parameters of the heat flow 
equations, C^, and R^, hive been discussed previously. 
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TUNNEL AERODYNAMICS 


The design of a closed circuit wind tunnel involves area-based, optimum 
circuit design of the endless tunnel duct so as to obtain desired flow 
velocity at the test section, assuring maximum flow quality at minimum power 
consmaption. The steady- state flow through the tunnel duct assumes different 
velocities around the circuit as dictated by the compressibility flow laws 
of continuity, momentum, and energy. The power losses in a closed circuit 
wind tunnel are caused by wall friction and vary as a function of cube of 
the speed; hence it is desirable to maintain the lowest possible velocities 
in the return ducting ocher than in the test section, which results in large 
areas of cross section and optimum aerodynamic streamlining. Circuit losses 
occur in various segments of the tunnel duct, viz fan nacelle, the forward 
leg diffuser, the third and fourth corners with turning vanes, the settling 
chamber, the screen section, the contraction, the test section, return leg 
diffuser, and the first and the second corner. Estimation of the circuit 
loss has been well studied (refs. 4 and 11). 

I 

The cryogenic tunnel is an endless pressure duct of varying cross- 
sectional aiea which has, in addition to regular components of a closed 
circuit tunnel, a liquid nitrogen injection section and a gas bleed section, 
located downstream of the test section and in the settling chamber area, 
respectively. Fan operation sets the gas in motion, and this flow in the 
endless duct can be a,nalyzed by forming the equations of flow continuity, 
momentum, and energy conservation for the tunnel. It is obvious that the 
flow process occurs both in space and time, a.nd appropriate distributed system 
equations need to be set up and solved (ref. 12) . 

The cryogenic tunnel is assumed to be made up of a large number of 
small length cylindrical segments of uniform area. Consider an ith 
cylindrical segment which is bound by sections i and i + 1 which are 
x(i) and x(i + 1) distant along the centerline from the reference fan exit 
section. Let the cross-sectional area of the segment be A(i) , mass flows 
entering and exiting be m(i) and m(i + 1) respectively. The total 
pressures at entry and exit be P(i) and P(i + 1), P(i) the average 
density, and T^(i) the average static temperature. The following discrete 


segments embracing the full 


equations set is valid for all i « 
length of the tunnel. 


Continuity 

The general equation of continuity for flow in a duct of cross sectional 
area A, fluid density p, mass flow rate m over a length 6x in a 
period 6t, is 

It 


This equation of continuity for the ith segment of the tunnel can be written 
in the discrete form as 


dt 


or 


A(i) p(i) [x(i + 1) - x(i)] = - m(i + 1) + m(i) 


Wg(i) = - ra(i + 1) + m(i) 


(17) 


In the cryogenic tunnel, under steady-state conditions, the mass flow difference 
m(i + 1) - m(i) is zero for all segments except for the ones having the liquid 
nitrogen flow inlet and gas bleed outlet. At these sections m(i + 1) - 
m(i) » m, or in respectively. Further, it may be noted that at steady state: 

U 5 


I 


^ W (i) =2^ m(i + 1) - m(i) 
i=l i*l 


' “g 


Momentum Equation 

The general equation of momentum associated with fluid flow in a cylin- 
drical segment of cross-sectional area A with a mass flow m, pressure P, 
velocity u over a length of 5x over a period 5t is 

A = It ™ - It 

This represents the force balance across the segment in the form of inertial 
forces of flow, the momentum flux induced forces and the wall friction drop 
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across. Utilicing the identify mfi) » p(i) A(i) uCi), the discrete form of 
the momentum equation can be written as 


A(i) PIL-t.-U..- P.(i) , . "(i) _ . m^Ci ^ D - m-(i) . 

x(i + 1) - x(i) A(i) p(i) A(i) p (i) [x(i + 1) 

Further, under steady-state operation of the tunnel 


x(i)] 


V P(i + 1) - P(i) * 0 


for the full tunnel circuit. If the fan segment is excluded 


^ P(i + 1) - P(i) = tunnel circuit loss AP 
i*l 

Energy Equation 

The general equations of energy associated with fluid flow in a cylindri- 
cal duct of length 5x and cross-sectional area A, fluid density p, static 
temperacure mass flow m, velocity u, and a metal wall that releases 

heat at rate of -Q^ can be written in a pair of equations 




^ cc ift T ) .. ^ 

5X P 3t 


-dg 


These equations correspond to the fact that total energy in the segment is the 
sum of internal energy and flow kinetic energy, and the energy rate equation 
includes temperature-dependent adiabatic heat rate, flow velocity contributions, 
and metal wall released heat. In the case of cryogenic tunnels, energy rate 
inputs also occur at the liquid injection section, the fan section and the gas 
bleed section. The discrete forms of energy equation can now be written as 






1 


^ {acD PCi) C^ T^(l) [xci*l) - x(i)] * - j - (i)^ii) 


■ Ts(i)]yCi) [X(i+1) - x(i)] ✓4irA(i) 


x(i+l? - 




[ra(i+l) 


ra(i)]T^(i) ♦ 


( 21 ) 


In the segments in which liquid injection, fan operation and gas bleed 

occur, additional energy rate terms of m.h, , Q- and -m C T are to be 

L L r 8 P 

appropriately added. 


FAN AND TUNNEL STEADY-STATE PERFORMANCE 


The steady state performance of a closed circuit tunnel and fan can 
now be analyzed (ref. 1), and values of mass flow, pressure ratio, fan 
temperature rise, and fan power at steady state determined. 


Mass Flow 


The test section mass flow in any tunnel is related to the settling 
chamber total pressure, total temperature and the flow velocity. Considering 
isentropic flow through contraction we have 


m * 0 u A » P„ AM 
s 


4 . 


Y T 


R mole 


Y + 1 


6965 ~ AM 

/r 


1 + 


Y-1 


m2 


2(Y-1) ]^g/ggc 


6965 — AM 

/f 


-3 


kg/sec , for y = 1-4 


1 + 0.2M2 


( 22 ) 


Pressure Ratio Across Fan 


In any closed circuit tunnel, the steady-state flow fan pressure ratio 
is a function of the pressure losses occurring in various segments of the 
tunnel and can be expressed in a simple form, normalized to test section Mach 
number as 


r = 1 + bM2 (23) 

where b is the integrated loss factor coefficient, and is likely to be a 
function of flow Reynolds number 


b = F(Re) 

with this knowledge of pressure ratio, the fan temperature rise can be 
estimated. 


Fan Temperature Rise 


The total pressure rise AP and the temperature rise AT are related 
under isentropic flow conditions as 


1 * ^ ^ , (i X lY-1 

P inlet ■* V T inlet/ 


Since the fan inlet temperature at steady state is T - AT, we have 


\ 


Y-1 


^ ' (t-at) ^ 


Simplifying 


AT * - 


Y-1 

rY - 1 


Y-1 

rY 


where n is the fan efficiency factor. 

Fan Power 

The fan power can now be estimated as the heat released by 


(24) 


• • 

^F ~ *^p ""fan inlet 


AT 


Though mass flow at fan inlet is m + ignoring the contribution of this 


small flow mj^. 


K_ P m 3 /f 

Qp r 3 

(1 + 0.2M2) 


(25) 


where 




^ AC b 
n P 


(^) 
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Spatial Profile of Tunnel Variables 


IVhen the cryogenic tunnel is in equilibrium, the steady-state solutions 
to the energy, continuity and momentum flow equations provide the spatial 
profile of the tunnel variables. Since c.(Uations (.I"’) and (20) are complex, 
the computations are simplified by describing the tunnel by 16 sections 
binding 15 segments as detailed in table 2. 


Mass flow: when the tunnel is in equilibrium 


then 


dW 

dt 


and m, = m 
L g 


m(9) » test section mass flow 

a 6965 — A(9) M (1 + 0.2 m2) 

/T 

m(i) * mC9) for 6 i ^ 10 

mCi) “ raC9) + m^^ for 1 £ i £ 5 

11 < i <_ 15 

Further, from the power point of view 

1 "^F P/f m3 

■"l "" ^ — * T3 

® B 3 (1 + 0.2M^) 

Tunnel Total Pressure 

The steady-state solution to the momentum equation is 

£ ^P(i + 1) - PCi)) = P(15) (r - 1) 
i=l 

The pressure profile can now be found to be 

P(9) = P 

i 

P(i] = P(9) - ^ 

9 


(26a) 


(26b) 
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The pressure drop AP in each segment is 


AP(i) . . AkCi) M^(i) P fi) (27) 

A^Ci) oCi) 

Values of estimated ^ ^ for the 0.3-m transonic cryogenic tunnel are shown 
in table 2. These values were obtained by distributing known tunnel circuit 
loss, based on simple analysis of geometrical and other loss effects. 

Total Temperature Profile 

IVhen the tunnel is in equilibrium, the energy rate of the system is zero. 
Based on this the tunnel total temperature profile can be determined based 
on equation (39) to be discussed later: 


T(9) 

» T 



i = 6 to 10 

T(i) 

= T - 

AT 

6+a 

6 

i = 11 to 15 

T(i) 

= T - 

AT 

C( 

0 

i ’• 1 to 5 


(28) 


Mach Number Profile 

The local flow Mach number M(i) at any segment i can be obtained by 
solving the sixth order polynomial in M(i) and seeking the solution between 
0 < M(i) < 1. 

[1 + 0.2M2(i)]^ - 6965 * 0 (29) 

m(i) /f(IT 

All quantities except M(i) are known in equation (29). 

Kinetic Energy 

An estimation of the flow kinetic energy associated with the gas flow 
through the tunnel can be made as follows: 


30 


( 30 ) 


E,(i) 


IVg(i) M2(i) YRTg(i) 

2j'g mole 


E^(i) 


, 338.9 PfU Vff, M^fi) 
mole • 2Jg 



Computer estimations of the various steady- state profiles and kinetic 
energy were made for various test section conditions of P, T and M, and 
these values are presented in the following text. 


Mass of the Gas in the Tunnel 

When the tunnel gas is stationary, the mass of gas in the tunnel is 

M, ■ 338.9 ♦ 250 — ] (31) 

8 T 1^ .^2j 

As the fan is accelerated, the density of the gas around the tunnel changes 
and density in the test section drops. As a result the total pressure goes up. 
Alternately, at a given pressure P, the mass of gas comes do\m. This change 
of mass has been function fitted as 

IVg = 338.9 ^ + 250 ^1 - 0.033 (32) 

Typical plots of spatial profiles of P(i) and Pg(i)> T(i) and M(i) for 
T = 100 K and P = 2 atm at test section Mach numbers of 0.3 and 1.0 are 
shown in figures 19 and 20. These profiles are plotted against tunnel center- 
line in the developed form. The strong correlation of cross-sectional area 
on the tunnel variables at different Mach numbers is evident. It can be 
further inferred that total and static pressures are grossly different only 
in the test section. Further, the flow velocity is very low in most of the 
circuit except at the test section. 

In figure 21, a plot of tunnel kinetic energy as a function of test 
section Mach number is presented. The kinetic energy is almost independent 
of the tunnel temperature and is of the type PM^ as detailed in equation 
(30) . At a tunnel pressure of 6 atm, the kinetic energy is maximum at 
about 190 kJ at M = 1. Comparison of the 0.3-m cryogenic tunnel kinetic 
energy with gas internal energy U shows that E^ is only about 

percent of U even at M = 1. Obviously, E^. is an insignificant fraction 
of metal enthalpy. 
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Tunnel Circuit Time 

The tunnel resident gas is moved around the duct in a finite time, which 
is the tunnel circuit time. Hence 

W * t 'in 
g c 

Ignoring the contributions of m^ entering the tunnel at section 11 and 
exiting at section 6, we have 


W 


^c ’ 


0.0045 


m(9) A(9) M»^ 


^ jl + 250 |l + 0.2M2 ^ j*! - O.OSSM^'Sj (33) 


An estimate of t^ was made for various P, T and M and is presented in 
figure 22. The tunnel circuit time is largest at low Mach number and low 
temperatures and smallest at Mach 1 and ambient temperature, 

Tunnel Fan Performance 

The fan performance can be expressed in terms of either the fan map or 
in terms of speed and pressure ratio data for the tunnel. Since the 0.3-m 
tunnel has been operational for a considerable time, a large quantity of fan 
speed, Mach number, and pressure ratio data is available. This steady-state 
data has been analyzed and function fitted to obtain following e.xplicit 
tunnel- fan interface performance. 

The fan speed and the test section Mach number are related by this 
emperical expression, and the fit is good to ilO rpm for any given Mach 
number : 


N = 597 M/f (1 - 0.3M) P'0.035 


(34) 


Secondly the fan pressure ratio and the test section Mach number are empiri- 
cally related as a function of flow Reynolds number. 


r = 1 m2 


[o.l97 (l - 2™ 


)] 


(35) 


The Reynolds number effects on r has been expressed in a very simple functional 
relationship. 
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TUNNEL TRANSIENT ANALYSIS 


The true transient response of the tunnel is the dynamical solution to 
the flow equations of continuity, energy and momentum both in space and time. 
In this analysis, a lumped model approach has been used in solving the flow 
equations by treating the tunnel as a single lumped volume. This can be 
justified by observing that flow kinetic energy which creates spatial vari- 
ation of tunnel energy states is a very small fraction of the total internal 
energy and the metal-stored enthalpy. These latter dominant energy terms 
dictate the tunnel temperature and pressure dynamics. However Mach number 
is a spatial function of geometry. Hence the approach used in the transient 
analysis has been to determine the average transient behavior at a low 
velocity segment on a lumped volume basis and superimpose the steady-state 
values of the tunnel spatial profile on this average state. The inaccuracies 
caused by the lumped approach are suppression of acoustic adiabatic waves 
and some circuit modes if the mixing is not good. 


Tunnel Dynamic Model 


Consider the lumped energy equation of a cryogenic tunnel system, from 
equations (19, 20) 


Surface 


A 

dt 


["gCvT, * * 




T 

g P 



Ignoring the rate contribution of kinetic energy and taking 
constant 


dW 


Surface 


“g s f * V \h - Vp"' -/ 


Surface 


If • rV - Cy”' * -/ '5t. 


to be a 


In order to allow evaporation and tran.sln time effects that are associated with 
each control input, three transport lag terms can be added. Since the tunnel 
conditions are measured at the settling chamber, transit times from the control 
input to the settling chamber are considered. However, the gas vent is located 
at the settling chamber and hence can be ignored; then 
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- 


Cv) T • q, 


-TpS 


Surface 

■/«t] 


( 36 ) 


The heat transfer from the metal wall occurs throughout the internal surface 
as was detailed in a previous section, This can be expressed as 


Surface ^ ^ ,j,g 

j () • y\ A 


where W^(i) is tie weight of metal in oach segment of the tunnel and t^(i) 
is the local metal time constant. Further, 


t (i) 
m^ ' 


where Yg(i) 


W^Ci) C, 


m 


[x(i + 1) - x(i)] /4ir A(i) 


2Yji) 


(37) 


1 


R„(i) 


» gas to wall heat transfer coefficient 


= 71.27 P(i)°*® T(i)‘0*22 M(i) [1 - 0.67 M(i)°-®5] 


Y (i) = n-T^ = (200 + 1.41 T) J/m2 K sec (ref. 13) 
m 

A computer evaluation of 't^(i) as a function of tunnel space has been made, 
and figure 23 provides the local metal time constant profile which indicates 
rapid response at the high-speed section and slow response at the settling 
chamber. Further, to lump the heat transfer effects, the following simpli- 
fication is made: 



C TS 
t m 

m 


where t is an average metal time constant for the tunnel 
m 


t = 
m 


w^c 

t m 


15 

E 

i*l 


[x(i + 1) - x(i)] /4n A(i) 


r2Y^Ci) Y^(i) 1 


2YgCi) . YJi) 


(38) 
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The gross metal time constant has been estimated for the whole tunnel envelope 
as a function of P, T and M and is shown in figure 24. Large metal 
time constants occur at low Mach numbers such as M ■ 0.2, P » l.S atm and 
T « 300 K. The lowest time constant occurs at M • 1, T ■ 100 K and 
P * 5 atm. Further, a simple function fit was made for use in simulation 
studies as 


t » 0.28 P“0*’ M’0‘7 sec 

m 


The tunnel temperature dynamics can now be e.xpressed as 
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The equivalent tunnel gas time constant t^ is a function of the relative gas 

thermal capacity to total thermal capacity and the metal time constant. This 

time constant t has been estimated as a function of P, T and M and is 
S JT 

show in figure 25. When = 0, equation 1591 provides the solution to 
spatial distribution of total temperature discussed in equation t2vSl. The 
total pressure dynamics of the cryogenic wind tunnel is contributed to by 
momentum equation and perfect gas properties 
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P - kaW T 


dt 


+ momentum effects 


A quasi-steady-state approach may now be taken, and the momentum effects are 
attributed to fan pressure ratio changes. In the fan 


fllL 

P(15) 
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1 + bM^ 


Assuming b to be a constant, Mach number changes and pressure ratio changes 
can be related as 


- - ■ (1 + bM2) + P(1S) [:b M 

Thus the total pressure dynamics is related to Mach number dynamics through a 
function of the type 


in 

dt 


= bMP 


dt 


For the cryogenic tunnel, the pressure dynamics can be expressed as 


dt 


P 

T 






+ K, . bMP 


dM 

dt 


(40) 


No transport delay terms are manifest in the pressure dynamics since pressure 
space dynamics occur at acoustic velocity. This transit time around the 
circuit is less than 50 msec at its largest. The Mach number at the test sec- 
tion can be directly related to fan rpm to provide the quasi-steady-state 
dynamics of Mach number 
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Because of the plenum volume, the Mach number dynamics is affected, and this 
change is of a first order type because of net flow into or out of the plenum. 
Then 

M « = N C41) 

k /f (1 + t S) 

m p 

where 

-V 

e = acoustic time delay 

tp = plenum time constant 

k = 597 (1 - 0.3M) P’0.035 
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LUMPED MULTIVARIABLE MODEL OF A CRYOGENIC TUNNEL 


The three equations (39), (40) and (41) derived previously provide 
explicit identities of average total pressure gradient P, total tempera- 
ture gradient f and the test section Mach number as a function of the three 
control inputs, the liquid nitrogen mass flow, the gas bleed mass flow and 
fan speed. The three equations represent the nine elements of a transfer 
function matrix model of a cryogenic tunnel which is shown in figure 26. 

This matrix model has linear dynamical terms with nonlinear coefficients. 
The parameters of this model are the tunnel geometrical constants, the 
thermophysical parameters of liquid and gaseous nitrogen, and the surface 
heat transfer parameters. Since the off-diagonal elements of the model 
are fairly dominant, the process is obviously strongly coupled. 

The nine input-output type gain elements of the multivariable lumped 
model of the cryogenic tunnel describe the steady-state process change for a 
unit incremental input of the control. The three diagonal elements repre- 
sent the direct gain terms from liquid nitrogen to temperature, fan rpm to 

Mach number, and gas bleed to pressure. The other off-diagonal terms are 

the interaction terms. The following study of the nine gain terms provides 
a knowledge of the dominance of inputs on the process and its coupling. 

In figure 27, a plot of temperature gain term (6 + aj/0 for a unit 
liquid nitrogen input is shown. Typically 1 kg of liquid nitrogen induces 

a steady-state change of gas temperature of -0.13 K/kg. However, under 

) 

transient conditions an approximate gain term of type (3 + a)/WgC^ describes 
the amplitude of response, which is shown in figure 28. The plot suggests 
a high temperature gain of -30 K/kg at P * 1 atm and T = 300 K, The 

gain reduces to -2 K/kg at P = 5 atm and T = 100 K. 

The term a/0 represents the temperature gain of the tunnel due to 
unit bleed of tunnel resident gas, and is shown in figure 29. This sensi- 
tivity is a small fraction of liquid nitrogen response and suggests a weaker 
coupling from bleed to temperature. The transient response is relatively 
higher than the steady state, but still weaker. 
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The tempei’ature gain term for fan speed is Qr,AN/N0 and is a very strong 

r 

gain element. Obviously the gain is high at high povcer, and the gain term can be 
approximated to ^ . At P * 5 atm, T * 275 K and M * 1 the gain 

is largest. 

The tunnel test section Mach number is directly related to fan speed 
and the coupling tei*ms from liquid nitrogen flow or gas flow are zero. How- 
ever the cooling terms from liquid nitrogen do affect the Mach number through 
the pressure ratio, which varies as a function of l/i-T. The Mach number gain 

for incremental fan speed is the constant k , which is nonlinear in Mach 

m 

number . 


The tunnel total pressure is influenced by all three control inputs and 
hence has corresponding gain terms. The liquid nitrogen affects the pressure 
by the terms t(P/T) and iP/hMAmj^. The gas bleed affects total pressure 
by the term (PW lAm , and the weakest influence is from fan speed through 

s s 

a gain term k. 


bPM 


df 


Thus, in the cryogenic tunnel model, seven of the nine gain elements are 
fairly dominant. Only the effect of gas bleed on .Mach number and the effect 
of gas bleed on temperature are low. Further, the direct effect of liquid 
charge on Mach number can also be ignored under regulated temperature con- 
ditions. 
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DESCRIPTION OF THE 0.3-m TRANSONIC CRYOGENIC 
TUNNEL SYSTEMS 


In this section, the engineering description of the 0.3-m cryogenic 
tunnel test facility and the various control subsystems relevant to modeling 
and control studies are provided (ref. 14) . 

0.3-m Tunnel 

The cryogenic tunnel is an endless duct made up of aluminum alloy 6061-T6, 
which is compatible for cryogenic operation. Tunnel area vs. centerline 
distance has already been detailed in table 1. The tunnel is made up of 
10 flanged segments which are coupled. The tunnel is anchored at the fan 
segment to a firm foundation, and the rest of the structure is mounted on 
counterbalanced supports and floating supports. The circumference of tunnel 
is 21.7 m with a minimum cross-sectional area of 0.1235 ra^ at tost section and 
a maximum cross-sectional area of 1.168 m^. The liquid nitrogen is introduced 
at a distance 15.85 m along the centerline from fan exit. The gaseous bleed 
occurs at a location 6.9 m from fan exit, the contraction cone is 2-m long 
and has a contraction ratio of about 10. The volume of the tunnel is 14.1 m^, 
and it weighs 3200 kg. The tunnel shell is externally insulated by four mats 
of chopped untreated fiberglass loosely sewn into mats of 7.6-cra thickness 
each and wrapped around with a fiberglass cloth which is purged with gas 
(fig. 30). 


Fan 

The tunnel shell is pierced by a fan shaft with an appropriate pressure 
seal and bearing, and the shaft is driven by a two-pole induction motor 
capable of delivering two MW of mechanical energy. This two-pole induction motor 
is driven by a rotary variable frequency generator with a frequency range of 
5 to 120 Hz and thus has a continuous speed control from about 300 to 7200 
rpm. The variable frequency generator can be remotely controlled from the 
control room. The fan shaft bearing is lubricated by a heated oil lubrica- 
tion system to prevent freezing of the lubricant. 

The fan motor system has been identified for its dynamics at optimum 
tuning, and the dynamics of the fan is 
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1 

I + o.sos + o.:s^ 


(4:) 


N 


N 


set 


The response of the system has near-critical damping and a natural frequency 
of 2.2 rad/~ec. 


Liquid .Nitrogen System 

Injection of liquid nitrogen requires a source pressure in excess of 
the tunnel pressure. In order to achieve this positive pressure, a liquid 
nitrogen pumping system is used. The schematic of liquid nitrogen pumping is 
shm^m in figure 31. The system consists of twin vacuum-insulated, stainless 
steel tanks constituting a total capacity of 212,000 liters. The storage 
system has a number of safety systems to maintain the boil off gases at a 
desired pressure of less than 3.5 atm (50 psiaj and backup lUowoff rupture 
disks to prevent overpressurication. 

The two tanks feed liquid nitrogen into a cryogenic, stainless steel, 
centrifugal pump driven by an AC motor. The pump capacity is 500 liter/min at 
10 atm pressure. Since the pump is basically a positive displacement device, 
the system pressure can be controlled by a pressure relief valve. This relief 
system is located nearly 50 m away at the tunnel. The 5-cm line carries the 
high-pressure fluid to the relief valve The relief valve is an air-to-open 
type pneumatic actuator driving a cry'Ogenic flow valve. The pneumatic actuator 
is driven from a three-term controller with a pneumatic set point generator. 

The line pressure is sensed by a pressure transducer and is converted into a 
pressure signal to be fed to the controller. This relief system can regulate 
the liquid pressure, and the pressure system dynamics has a 6 rad/sec band- 
width and O.S damping. The output of the relief valve is returned to the 
top of the tank through a 50 m long 5 cm line. The tunnel is fed from the 
regulated section of the liquid nitrogen line through four digital valves. 

Liquid Nitrogen Control Valves 

The liquid nitrogen flow into the tunnel is controlled by four digital 
valves. The digital valve has a number of solenoid-operated elements which 
are pure binary weighted. The valve flow coefficient is 4. The smallest 
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element allows 1/1024 of full open> and subsequent elements progressively 
double their area. All the elements are solenoid operated and have a sonic 
nozzle. The flow equation is 


m 


L * 56.805 \ 


or 


iTij^ a 5.47 - P kg/sec 


0 < Aj^ < 1 

step size * 1/1024 

3-47 ^ (** 3 ) 


for a liquid nitrogen density of 0.804. 


Gas-Bleed Valves 

The cryogenic tunnel pressure shell is certified for six atm. To prevent 
the possibility of overpressurization, the settling chamber has three large 
bleed valves open to atmosphere, each controlled pneumatically. The valves 
are air-to-close type to allow fail safe shut down. Each valve can be sepa- 
rately controlled from the control room. One of these valves is cascaded 
with a 10-element digital valve meant for closed-loop operation. In view of 
the large areas necessary to control gas flow at low pressure differentials, 
the elements of the digital valve are air-piloted solenoid valves. The valve 
flow areas are binary weighted sonic nozzles. 

Such a valve has a flow, which under choked conditions, is 


""g “ V If ^^*^d 


Assuming a C* = 0.7 and a valve coefficient of 8, 


m * 21.8 -i- A„ kg/sec 
g /f g 


0 ^Ag £ 1, step size = 1/256 (44a) 


when the pressure ratio across valve is very low unchoking occurs. Then 


m = 21.8 A„ 
g /f g 


■> ■ (¥)'"■ 


kg/sec = k„ — A„ kg/sec 
g /f g 


for 1 < P < 1.5 


(44b) 
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Summary of Model injj Effort 


Hitherto the svTithesis of a lumped multivariable model of a closed- 
cii'cuit, transonic, cryogenic pressure tunnel with specific reference to the 
0.3-m cryogenic tunnel was discussed. The three input/three output model 
so generated describes the tunnel variables in terms of tunnel physical 
parameters, thermophysical properties of nitrogen, and the tunnel control 
input in an explicit manner. The tunnel model embraces the global envelope 
of operation and demonstrates the interactive nature of the process and its 
nonlinearities . 

In deriving the lumped multivariable model of the cryogenic tunnel, 
obviously some of the faster d\'namic modes of the tunnel have been suppressed. 
A heuristic study of the cryogenic tunnel flow phenomina indicates possible 
dynamic modes described by eigenvalues attributable to turbulance modes, 
acoustic adiabatic modes, two-phase flow induced modes, circuit modes, gas 
and metal thermal modes, and the slow mass-energy integration modes, in a 
descending order in frequency. In the lumped model svTithesiced here, the 
faster modes like circuit, acoustic, and turbulance modes have been ignored, 
because they constitute the uncontrollable modes for the degrees of freedom 
available at conti'ol inputs. 
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VALIDATION OF THE MODEL 


A model of a physical process can be validated only by comparing it with 
the actual process. Such a comparison can be at parametric level or at a 
performance level. In the former, the actual process records can be analyzed 

using modern identification procedures to arrive at tunnel gain parameters 
and transient parameters. The parametric identification of process is an 
involved procedure which requires input disturbance design and experimental 
record analysis using identification procedures. Identification of nonlinear 
systems is difficult. In the case of comparison of performance, the model 
can be used to predict the process response for various types of inputs, and 
the actual tunnel responses and the predicted responses can be compared for 
quasi-steady-state and transient responses. 

In the present research effort, the latter route of comparing time 
responses of the tunnel to the predicted time responses of the model was 
chosen. In order to generate the model responses, a real-time simulator was 
developed which is a true electrical analog of the process equations and 
model derived in the previous section. In the following sections, the details 
of the development of such a real-time cryogenic simulator are discussed. 


SIMULATION’ OF THE 0.3-m TRANSONIC CRYOGENIC TUNNEL PROCESS 


The methods available for simulating any process in real time are on 
analog, digital or hybrid computers. Analog computation is very fast, but 
its accuracy and its ability to perform accurate nonlinear computations are 
very limited. Digital computation is very accurate and its speed because of 
serial operation is limited by the quantity of computation. Hybrid analog- 
digital computers with dedicated digital processors provide an efficient means 
for combining speed and interactive features of the analog computer with the 
accuracy of digital computers. In the present effort, the simulator of the 
0.3-m cryogenic tunnel process has been carried out on a hybrid computer. 

Hybrid Computer 

The hybrid computer system used to develop the 0.3-m transonic cryogenic 
tunnel simulator consists of (1) a 16-bit, 32-k memory digital processor with 
a cycle time of 1 usee and a floating point hardware for arithmetic, (2) a 
10-V, 100 amplifier analog processor with logic and nonlinear units, (3) an 
interface hardware unit consisting of analog to digital converters and 
digitally controlled attenuators, and (4) an interface system between the 
processors for control, logic, and high-speed data interface. 

The hybrid computer has I/O units in the form of a card reader, line 
printer, a CRT control unit, a cassette tape unit, and an analog recording 
system. The real-time operating system is through a real-time operating 
system and allows use of real-time FORTRAN and assembler. 

Simulator Display and Conlrol Unit 

In order to achieve real-time interactive simulation of the 0.3-m cryo- 
genic tunnel to allow man- in- the- loop operation, a control display unit 
similar to the actual tunnel control was made. The operator controls consist 
of a fan speed control, liquid nitrogen valve area control, and gaseous nitrogen 
bleed area control. The operator displays are gas and average metal tempera- 
ture in Kelvin, test section Mach number to 0.001 resolution, and tunnel total 
pressure in psia (0.07 atm steps) with a resolution of 0.1 psia (0.00" atm). 

The operator display and control unit is shown in figure 32. The other displays 
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on this unit are the fan speed in rpra and the gas valve and liquid valve 
areas in percent of full opening. 

The hybrid scheme used for simulation has been to compute nonlinear terms 
on the digital processor and the linear dynamical terms on the analog proces- 
sor. The digital computer work load is shown in figure 33 along with the 
hybrid scheme. The IS listed computations are the typical ones and are car- 
ried out in about 35 msec. 

The numerical inputs to these computations are derived from the outputs 
of the ADCS appropriately scaled into engineering units compatible with the 
process equations (36) through (44) . The computed outputs are rescaled into 
analog processor units and are fed to the analog processor through digitally 
controlled attenuators. This cycle of analog and digital computation update 
occurs at a rate of 2S Hz. 

The analog processor handles the linear dynamic terms in real-time scale 
and amplitude scaling based on static quantities. Only the temperature 
derivative f is rescaled in its amplitude. The analog processor runs contin- 
uously in real time and its coefficients are updated by the digital processor 
once every 40 msec. The appendix provides the typical digital processor program 
used. Figures 34 through 37 provide the detailed analog processor patching. 

The temperature dynamics of the tunnel is patched as shown in figure 34. 

The gas temperature is obtained by integrating f, which in turn is obtained 
from subtracting Q^/WgC^ from f + Q^/W^C^. The temperature data is boot- 
strapped to generate metal temperature T^ using time constant t^. Figure 
34 details generation of gas temperature, metal temperature, and the thermo- 
couple output. 

The analog scheme for generating pressure, fan and Mach dynamics is 
shown in figure 35. The total pressure P is derived from integrating P 
and accomodating the fan speed coupling. The second order fan dynamics is 
shown and uses two integrations with the N integrator being limited to -600 
rpm/sec acceleration. The test section Mach number dynamics is a simple first 
order filter. 

In figure 36, the details of the liquid and gas valve drive are described 
and the provisions for manual control, transient pulsing with known amplitude 
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and period, and closed- loop control are shown. In figure 37, the closed-loop 
controllers for pressure and temperature loop are shown. These loops are 
the proportional-integral type and have loop gain scheduling built into them. 

Thus figures 33 to 37 detail the analog segment of the hybrid computation. 

The analog segment has provisions for creating initial conditions for P, T 
and M; The analog processor has continuous performance, and 10 analog quantities, 
vie P, T, M, N, A^, Ag, 30 f, t and T^, are AD-converted and 

fed to the digital processor. Thus, totally the system represents a hybrid 
scheme for simulating the 0.3-m cryogenic tunnel. 

The digital processor program was written in REAL F0RTRA.N and assembler 
(. Appendix! , and the system was run continuously, updating being at 25 He. 

The bandwidth of the process is no more than 1 to 2 Hz, and hence the system 
is truly continuous. 

The system response were recorded on a 40 -sample/ sec digital recording 
system for processing and further analysis. 


» 
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PERFORMANCE TESTS 


Three types of validation performance tests were performed, both on the 
cryogenic tunnel proper and its simulator: steady-state equilibrium tests, 

quasi-steady Iona-term warmup/cool down tests, and transient response tests. 
The tunnel and the simulator performances were then compared to check for 
validation of time responses. Further, the responses were correlated with 
the math model terms. 


Transient Response Tests 

It is well known that linear time- invariant systems yield characteristic 
responses to deterministic input disturbances. The cryogenic tunnel is time 
invariant but has a nonlinear, complex structure. However, by assuming local 
linearity the system can be made to yield transient responses which can be 
interpreted on the basis of linear system responses. The cryogenic tunnel 
model has integration type dynamics in pressure and temperature. Hence step 
inputs cannot be used because of resulting divergent responses. Sinusoidal 
inpits are difficult to realize in respect to liquid and gas flow. Hence 
impulse response testing has been resorted to. 

During the inpulse response tests, both the 0.3-m transonic cryogenic 
tunnel and its simulator were first balanced at specified equilibrium condi- 
tions by open-loop adjustment of the three coiitrol inputs. Then the tunnel 
and the simulator were subjected to impulse disturbances of liquid nitrogen 
flow, gaseous nitrogen bleed flow, and fan speed at successive intervals. In 
order to study the dynamic behavior of the tunnel and its simulator over the 
entire range of tunnel operational envelope, the tunnel variables P, T and 
M were scanned from P = 1.5 to 5 atm, T = 100 K to 275 K and M = 0.3 to 
0.75/0.9 in 27 different combinations as in table 3. 

The amplitude and the period of the pulse were chosen so as to produce 
easily measurable changes in the tunnel variables. Thus a 6.25 percent of 
full area increment to the liquid nitrogen and gas bleed valves and a 100 rpm 
decrement to fan speed were chosen for the pulse amplitude. The period of 
pulse was uniformly 3 sec, which corresponds to 0.75 to 4 circuit periods. 
Typical responses to these pulses were 0.5 to 6 K in T, 0.016 to 0.1 atm in P 
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and 0.006 to 0.03 in M. A 30-s®c settling time was allowed between pulses. 
Tests on the tunnel and the simulator were performed separately. 

The tunnel and the simulator responses for identical tunnel conditions 
are presented in matched pairs in figure 38 to 64. These correspond to 
tunnel conditions in table 3. The tunnel responses were generated by a manual 
operator induced 3-sec pulse whose period was only approximate. The liquid 
nitrogen pressure also fluctuated between 10 to 20 percent off the nominal. 
This prevented precise knowledge of either Am, or Am , the incremental 

M g 

flow. An inspection of the figures 3S to ‘ indicates noise-free responses on 
the simulator and noisy responses on the tv ael, particularly in respect to 
Mach number. 

The typical impulse response of the tunnel is presented as a response 
trajectory A through K for the total pres.sure test section Mach number and 
total temperature. These response trajectories A through K arc now analysed 
for each control input pulse. 

Liquid Nitrogen Pulse Input 

The trajectories ABCD, of tv'pical responses in figure 38 (55 or So), 
correspond to tunnel dynamics in P, T and M for a liquid nitrogen pulse 
input. Prior to point A, the cryogenic tunnel is in equilibrium with 
^ = f = 0 and m, « m . During period AC of the trajectory, an extra liquid 
flow Araj^ is imposed on the tunnel. 

The tunnel gas temperature T reduces in period AC, due to negative T 
as in equation 39. Since the wall temperature lags the tunnel gas, T is 
largest at A and progressively reduces in period AC as dictated by the lead 
lag terra. This effect is due to wall heat release described by metal time 

constant t . At point C. the extra liquid nitrogen is cut off. The gas 

S , 

temperature gradient T is now positive because of metal being warmer than 
gas. Hence gas temperature T increases in period CD, and at point D bulk 
of heat transfer effects are complete. The metal and gas temperatures sta- 
bilise. A transport delay between gas temperature and command, corresponding 
to particle transit time from liquid injection station to settling chamber, 
is noticeable. 
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The gas temperature at point C is lowest and AT away” from the starting 
value. This AT can be related to transient temperature gain of figure 28. 
Study of tunnel and simulator responses to liquid nitrogen pulse, of figures 
38 to 64, generally confirm the relations of figure 24. At low gas weight 
(figs. 56 to 581 AT is very large, and at high gas weight (figs. 44 to 46) 

AT is smallest. However discr*epancies between tunnel and simulation responses 
are evident in figures 56 to 58. Since the tunnel responses are fi'om a thermo- 
couple, the differences can be attributed to first order filtering that occurs 
in this sensor. 

The total pressure trajectory ABCD can now be analysed (figs. 38, 55, 
561. .4s detailed in equation (40), the total pressure gradient P is the sum 

of a cooling terra ^ ® '"*^ss term Am^^. During period AB the cool- 

ing term is greater than the mass term and® P is negative. .At point B, the 
two are equal and pressure is at a minimum. During period BC the mass terra is 
greater than the cooling term and P is positive despite f being negative. 

At point C the excess liquid flow is made zero and P continues to be positive 

P • 

because of term y T. At point D pressure stabilizes on completion of bulk 
of heat transfer. 

Thus a typical pressure response of the cryogenic tunnel to a liquid flow 
pulse is characterized by pressure falling initially and after reaching a 
minimum the pressure starts rising during the pulse period. All 27 pressure 
responses, of figures 58 to 64, reveal this signature of the tunnel. 

In case of low temperature 100 K responses, of figures 33 to 39, the 
minimum pressure point consistently occurs at or very near point C in both 

T 

the tunnel and the simulator. Thus at low temperatures cooling effect P t 

T 

dominates the mass effect because of relatively poorer heat release from the 
walls. It may be noted that at 100 K metal enthalpy is 12 percent of that at 
300 K, as in equation (11). Hence pressure rises over the period BD are 
relatively small. 

In contrast, the high temperature runs at 275 K (figs. So to o4) consist- 
ently locate the minimum pressure point B at or very near A, indicating a 
P • T 

dominant Am. over the cooling term P — . Hence the tunnel pressure 
♦V u 

g T 

rises almost from the instant of liquid injection. The pressure rise in period 
BD is comparatively larger because of higher met.al enthalpy, lower gas weight, 
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and high specific heat of metal [eqs. (9) and (10)]. The tunnel and simulator 
responses demonstrate this feature. 

The 200 K runs (figs. 47 to 55) show a minimum point B occuring at all 
phases of AC. The point B shifts from near A (as in fig. 44) to C (as in fig. 
48) with increasing gas weight and reduced Mach number (corresponding to lower 
convection heat transfer). 

The pressure gradient at the peak pressure near point D in the post-LN 2 

pulse period makes an interesting point of study. At its peak near point D, 

the gas bleeding is greater than m^ [eq. (40)] and the pressure gradient is 

a function of A . The steady state A [eq. (26b)] exponentially increases 
8 ^ S 

with Mach number. Hence largest P occurs at M = 0.75/0.9 as in figures 

49, 52, 55, 58, 61 and 64, and P is smallest at M = 0.5 as in figures 47, 

50, 53, 36, 59 and o2. However, since peak pressure P is very nearly the 
same as starting at 100 K, these effects are not very evident in 100 K runs 
(figs. 38 to 46). Both the tunnel and the simulator confirm the features 
discussed. No transport delay is evident in the response. 

f 

The Mach number trajectories ABCD of figures 38, 55 and 56 can now be 
analyzed. In the period AC, the Mach number rises linearly and reaches a 
maximum at C. This can be explained as being due to change in the fan pressure 
ratio caused by cooler gases passing through the fan. The fan performance 
at constant speed can be noted to be of type M/f and is constant [eq. (41)]. 
Since the tunnel gas temperature drops in AC, the Mach number rises. Simi- 
larly, during period CD the extra liquid is made zero and temperature rises; 
hence Mach number drops linearly to the original value. 

At point C, the Mach number is largest and has an incremental value 
AM. This response amplitude can be related to gas temperature by equation 
- 'Iy derived from equation (41). The actual AM in all the runs of 
figures 38 to 64 can be used to estimate AT. Typically, in figure 38, the 
test section Mach number is 0.30 with AM = 0.006 at 100 K, which corresponds 
to AT of -4 K. However, the recorded AT is only -1.5 K, and the 
difference can be attributed to attenuation caused by the first order dynamics 
of thermocouple. In all the 27 cases, the tunnel and simulator responses 
show good agreement. No transport delay occurs since the coupling is through 
fan. 
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In summary^ the response trajectories of the tunnel variables to a liquid 
nitrogen pulse correlate well with the math model responses predicted by the 
simulator. Pressure responses with a minimum occuring during pulse period 
show a characteristic signature of a cryogenic tunnel with metal in the tunnel. 
The model and tunnel responses agree quite well over the entire range of tun- 
nel operation. The test section Mach number responses also show generally 
good agreement throughout. Relatively poorer agreement in gas temperature 
responses are evident in some cases. This can be attributed to thermocouple 
dynamics. Since Mach number match between the tunnel and the model are good 
and temperature predictions based on Mach numbers do agree well, it appears 
that the thermocouple response varies as a function of P, T and M. 

Further, transport delay effects are evident in all the tunnel and model 
responses and agree generally. 

Gaseous Nitrogen Bleed Pulse 

The response trajectories EFG correspond to the tunnel dynamics in P, 

T and M for a gas bleed pulse input. Prior to point E, the cryogenic tun- 

« • • • 

nel is almost in equilibrium with P = T = 0 and m = m. . During period EF, 
the gas valve is incrementally opened to allow an extra flow Am out of the 
tunnel . 

The tunnel total pressure P drops linearly with a negative P. This 

P . p • 

gradient is caused by mass term ~ Am^ and by a second order term. The 

contribution of the cooling term if very small. The pressure P reduces 
linearly in the period EF. At F the incremental valve area is made zero 
and Aiig » 0. However, the pressure at point F is lower than the starting 
valve and hence m^^ is greater than m^, resulting in a positive gradient 
for P during FG. 

The pressure at F is lower than the starting pressure by a magnitude 
AP, which is a function of AQ and in turn is proportional to Pv^. 

The largest gradient occurs at P = 5 atm and T = 275 K as in figures 62 
to 64, and the smallest gradient occurs at P = 1.5 atm and T - 100 K as 
in figures 38 to 40. The tunnel and the simulator responses agree well 
throughout as in figures 38 to 64. 
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During period FG, the pressure gradient is positive and is a function 
of steady-state valve opening Ag. However the valve Ag is a function of 
fan power and is largest at M ■ 0.75/0.9. Thus the pressure recovery is 
fastest at M ■ 0.75/0.9 as illustrated in figures 40, 43, 46, 49, 52, 55, 
58, 61 and 64 and slowest at M ■ 0.3 as in figures 38, 41, 44, 47, 50, 

53, 56, 59 and 62. No transport delays are evident in the tunnel pressure 
responses . 

9 

The tunnel gas temperature T shows a very poor coupling to gas bleed. 
Energy equation (39) predicts a negative T for gas bleed. This gain 
sensitivity is shoivn in figure 29 and is comparitively very low because 
a < 6. During the period EF the temperature tends to drop but is filtered 
by the lead lag term due to wall heat release. The temperature T at F 
can be appro.ximated to AQ . Since is a function of P//f, AT 

"'g^V ® ® 

is only a function of temperature T, the largest responses occur at 275 K 
and smallest at 100 K. Figures 38 to 64 generally confirm this, but many of 
the responses are small in amplitude. The temperature recovers slowly in 
period FG to its original value because of continued he.it release. 

The tunnel Mach number responses to gas bleed pulse are insignificant. 
The coupling between Mach number and gas bleed occurs dominantly thi'ough 
temperature only. The test section Mach number is related to pressure by a 
very weak empirical function p0.035 be ignored [eq. (41)] for small 

perturbation. In principle, during period EF the Mach number ought to rise, 
demonstrate a peak at F, and fall back to original value by B. In the tunnel 
responses, the noise submerges these changes. In the model responses which 
are noise free this effect can be observed in figures 38 to 64. 

In summary, the tunnel responses to a gas-bleed pulse correlate well 
with the model predictions in the form of simulator responses. The most 
dominant effect of gas bleed is in respect to total pressure, which falls 
linearly as predicted. The temperature responses to the gas bleed pulse are 
quite small, and the simulator predictions agree generally. The Mach 
number responses of the tunnel to gas-bleed pulse are very insignificant. 
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Fan Speed Pulse Input 


The response trajectories HJK correspond to the tunnel dynamics in 
P, T, and M for a fan speed pulse input. Prior to point H, the cryogenic 
tunnel is almost in equilibrium with P » f » 0 and liij^ » m^. 

During period HJ, the fan speed command is reduced by 100 rpm and the 
fan decelerates as defined in equation (42) in a second order fashion. The 
decreasing speed results in reduced fan pressure ratio, which in turn mani- 
fests itself as decreasing Mach number profile HJ. During the period JK, 
the fan is accelerated back to original speed by restoring full speed. The 
fan pressure ratio increases gradually in JK as does the Mach number and is 
restored to near original value at K. The Mach number response of the 
cryogenic tunnel to a fan speed pulse is almost a linear fall in HJ and a 
linear rise in JK with maximum deviation AM occuring at J. This AM 
is related to AN//T as in equation (41) . Since AN is constant for all 
the runs, AM is proportional to l//f The largest decrement of about 

0.034 occurs at 100 K (figs. 38 to 46), and the smallest occurs at 275 K 

(figs. 56 to 64). The model and the tunnel Mach number responses match rea- 

sonably well throughout. The Mach number response to fan speed pulse has 
another constituent via temperature drop which opposes the pressure ratio 
effect. However, its contribution is very small. The Mach number responses 
to fan speed do not indicate any transport delays. 

The tunnel total temperature trajectory HJK, due to fan speed drop, 
shows a reducing temperature T in the period HJ and a very slow rise in 
JK. This is due to decreasing Qp in HJ, which results in a negative 

T, which is affected by wall heat release, as in equation (37). The 

maximum deviation occurs at J, when the original fan speed command is 
restored, and Qp reaches normal value at K. During period JK the gradient 
f is changed from negative to positive due to wall heat release. 

The fan power drop AQp is maximum at J and is proportional to 
Pv^T f (M) as in equation (25) . Thus maximum AT occurs when AQp is 
maximum as at M = 0.75/0.9, P = 5 atm, and T = 275 K as in figure 64. 

The lowest response occurs at M = 0.3, P = 1.5 atm, and T = 100 as in 
figure 38. The trajectories HJK in temperature show a time delay from 
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command to response corresponding to particle transit time from fan to set- 
tling chamber, which is a function of tunnel circuit time. The tunnel and 
model responses agree with each other in their temperature responses gener- 
ally (figures .38 to 64) and confirm their relation to AQp. 

The tunnel total pressure trajectories HJK, due to fan speed pulse, 

can be analyzed using equation (40) . The pressure gradient is contributed 

P • dM * 

to by =r T and by k^., bMP gr* During the period HJ, T is negative and 

dM . , • • 

-j— IS negative. Hence P is negative and pressure continues to fall during 
HJ. At J the fan speed command is restored and is positive with f 

slowly reaching a positive value in the period JK. The maximum pressure 
deviation AP occurs at J. Since the effect ^ is most dominant, AP 
can be related to product MP. Thus largest AP occurs at P * S atm and 
M » 0.75/0.9 as in figures 46, 55 and 64. The smallest AP occurs at 

P a 1.5 atm and M = 0.3 as in figures 38, 47, and 36. The pressure re- 

sponses do not show any transport delay features to fan speed command. 

In summary, the fan speed pulse on a cryogenic tunnel creates a dynamic 
change in all tunnel variables M, P and T. The Mach responses are most 
dominant and are functions of 1//T. The pressure and temperature responses 
are functions of fan power and are evident only at high P and M. Transport 
delay effects are observable only in respect to tunnel temperature.. The 
tunnel responses generally agree with the model responses predicted by the 
simulator. 

Tunnel responses of figures 39, 52 and 57 need clarification. Oscillatory 
responses of figure 39 were caused by inadvertent cycling of one of the 
smaller elements of the gas valve. The fan rpm command in figure 52 actually 
consisted of a 0.5-sec pulse followed by a 3-sec gap and another 3-sec pulse. 
All the tunnel variables show the cumulative effect of twin pulses. In figure 
57, the large temperature response of the tunnel for LVo flow was caused by a 
five-sec LN 2 pulse. 


QUASI -STEADY-STATE TESTS 


In addition to the transient response test, studies on the cool down 
of the tunnel from ambient and warmup of the tunnel to ambient, using excess 
liquid nitrogen flow and excess fan heating respectively, were conducted in 
order to validate the long-term dynamics of the heat transfer model. 

The basic energy equation of the cryogenic tunnel process is 


Surface 

dT i* *v* • • /•• 

Vv 3t * (™L ■ ""g) * ■ Vp’^ -/ 


where 


Surface W C T 

-\)y 


m 


In order to prove this heat transfer model, time histories of average metal 
temperature T^^ and gas temperature T were recorded for ordered energy 
input terms. Matches between temperature, metal to gas temperature difference, 
and time were sought for the tunnel and simulator quasi-steady-state responses. 


Cool Down Studies 

The tunnel cool down test consisted of starting with a metal and gas 
temperature at ambient of 290 K and running the tunnel with excess liquid 
nitrogen flow while the tunnel pressure and low fan speed were maintained 
constant . 

Under these conditions m, “ m because of constant pressure and 

^ o 

Qp « ' ^p^) speed, then 

Surface 

Yv ^ 

«t ■ {'f - \) y ■ ’’n, ‘ TTTT 

m 

The tunnel and the simulator were balanced at an equilibrium condition of 
290 K, 1.5 atm at 1200 rpm. The liquid nitrogen pressure was maintained at 
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t) atm. The tunnel cooling was started by rapidly opening the liquid valve 
to 12.5 percent open, corresponding to about 1 kg/sec flow. The tunnel 
pressure was maintained at 1.5 atm throughout. Gradually the tunnel metal 
cooled do\m, and time histories of gas temperature, me cal temperature at a 
forward diffuser leg location Cfig- -3) and the test section Mach number 
were recorded. Figure 65 shows the time histories of tunnel gas and average 
metal temperature, the tunnel gas and metal pressure. The tunnel cooling 
was taken do\vn to about 100 K in about 35 min. It may be noted that model 
and the tunnel temperature histories match reasonably well. An initial gas 
to metal temperature difference of about 65 K (at 290 K) and a final dif- 
ference of 15 K (at 100 K) occur in both the tunnel and simulator. The simu- 
" lator and tunnel tempei'atures agree to within 5 K over most of the cooling 
period e.xcept at the very end where the difference is 10 K. This can be at- 
tributed to heat gains thi'ough some flanges in the actual tunnel which make 
cool dowTi a bit slower. The Mach number time histories of tunnel and simulator 
agreed to better than 0.005, starting from an initial value of 0.140 and ending 
with a value of 0.228. The total quantity of liquid nitrogen consumed for 
cool down from 290 K to 105 K was determined for both the tunnel and simulator. 
The tunnel value was 1825 kg against a simulated consumption of 1800 kg. 

Thus, a fairly good match between the tunnel and the model for a quasi-steady- 
state cool down confirmed the heat transfer modeling used. 

Warmup Studies 

The tunnel warmup test consisted of starting with a cold tunnel with 
both the gas and metal at 100 K, and then running the tunnel at fairly high 
Mach number with the liquid flow fully cut off and the tunnel pressure 
maintained constant, then 


m, “ m =0 
L g 


and Q„ >> m C T 
F g P 


-IK 


Surface 


g 


. . W. TS 
In t m 

^t ’ 1 + t S 

m 

Kp PM3^T 

(l + 0.2 
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Both the tunnel and the simulator were cooled to 100 K and were balanced at 
100 K, 2 atm, and a Mach number of 0.6. 

The warmup was started by cutting off the liquid nitrogen flow into the 
tunnel, then mainc...ining the tunnel pressure and Mach number by fine adjust- 
ment of the gas valve area and fan speed. Records of gas temperature and the 
average metal temperature were taken as functions of time. The average metal 
temperature in the tunnel occurs and is sensed at a location in the forward 
leg diffuser (see fig. 23) , which provides a metal time constant as a function 
of tunnel space. The warmup profile of the tunnel and the simulator at 
M » 0.6, P * 2 atm from 100 K to 300 K is shown in figure 66. The metal to 
gas temperature differences agree very well. However, the tunnel shows an 
initial warmup rate lower than the simulator, and this has been attributed to 
an initial tunnel pressure of less than 2 atm. Subsequently the tunnel and 
simulator warmup gradients match quite well. The metal to gas temperature 
difference in the tunnel and the simulator agree to within 2 K, indicating 
the adequacy of the heat transfer model. 

Validation 

In summary, the 0.3-m transonic cryogenic tunnel and its model responses 
show acceptable agreement, not only in their transient responses to various 
control inputs over the entire range of operational envelope, but also in 
their long-term quasi-steady-state warmup and cool down characteristics. All 
the features of the transient and quasi-steady responses correlate well with 
math model predicted effects. 

This overall steady-state, transient, and quasi-steady long-term re- 
sponse match leads to the conclusion that the model is a basic representation 
of the tunnel and that the model is validated. 


CON'CLUDIXG REMARKS 


In this report, the development of a conti'ol compatible mathematical 
model of a liquid nitrogen cooled, gaseous nitrogen operated, transonic, 
cryogenic pressure tunnel has been detailed. Utiliting the thermophysical 
properties of liquid nitrogen, gaseous nitrogen, and the tunnel metal 
aluminum, the thermodynamics and the flow dynamics of the cryogenic 
tunnel have been analyzed. The concept of controlling the energy state 
of tunnel gas provided a simple understanding of the energy associated with 
various physical processes occuring in the tunnel. The metal enthalpy was 
proven to be the most dominant store of energy in the system followed by 
the internal energy of the gas. The flow-induced energy terms, otiier than 
wall friction, turned out to be a very insignificant fraction of the total 
energy of the system, thus providing an idea of i*elative importance to be 
attached to various energy storage terms. 

The basic control inputs to the tunnel were identified as rate vectors 
of gas mass and energy either entering or leaving the system. The tunnel 
insulation makes the system an autonomous thermodynamic system with mass 
and energy crossing the boundary. This open thermodynamical system was ana- 
lyzed in order to arrive at a very simple, lumped, multivariable, global, 
mathematical model of any closed circuit, cryogenic pressure tunnel. The 
parameters of this model are nonlinear functions of thermophysical parameters 
of the tunnel metal and gas and tunnel variables. 

Secondly, validation of the lumped, multivariable model using the 
tunnel transient and quasi-steady-state responses was discussed. The 
validation has been based on matching the impulse responses of the tunnel 
with the model. The model responses for various disturbances were derived 
using a cryogenic tunnel simulator. The cryogenic tunnel simulator is an 
interactive real-time simulation of the various model equations realized 
on a hybrid computer. The hybrid computer uses the digital processor to 
estimate comple.x nonlinear product terms, which updates a continuous running 
analog computer. 

The validation has been demonstrated by good agreement of both transient 
and quasi-steady responses between the tunnel and the simulator under 
similar operating conditions. 
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Table 1. 0.3-m TCT area distribution. 


DISTANCE, X 

AREA, A(x) 

DISTANCE, 

. X 

AREA, A(x) 

(m) 

(m2) 

Cm) 


(m2) 

0 

0.254 

12.440 


0.267 

0.635 

0.282 

12.540 


0.2192 

1.245 

0.357 

12.640 


0.184 

1.854 

0.437 

12.750 


0.158 

2.641 

0.516 

12.85 


0.138 

3.251 

0.566 

12.950 


0.129 

4.470 

0.816 

13.050 


0.1250 

5.690 

1 . 006 

13.220 


0.1235 



- 

TEST 

SECTION 

6.299 

1.108 

14.14 


0.1250 

6.908 •<- GNo 

flow 1.168 

14.25 


0.130 

8.483 

1.168 

14.52 


0.145 

9.956 

1.168 

14.65 


0.152 

11.328 

1.162 

15.19 


0.189 

11.430 

1.149 

15.64 


0.210 

11.530 

1.112 

■ 15.850 

^ LN2 

flow -> 0.240 

11 . 630 

1.042 

16.307 


0.284 

11.730 

0.948 

16.764 


0.330 

11.830 

0. 833 

17.221 


0.379 

11.930 

0.719 

17.678 


0.4338 

12.030 

0.597 

18.135 


0.4552 

12.140 

0.496 

19.354 


0.455 

12.240 

0.404 

20.570 


0.394 

12.340 

0.328 

21.285 


0.253 




|Fan| 




21.540 


0.252 
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Table 2. 0.3-m TCT circuit loss pattern. 


DISTANCE 

(m) 

0.00 

0.521 

1.435 

2.458 

4,779 

7.214 

8.759 

11.061 

13.075 

14.630 

16.531 

15.105 

19.627 

20.080 

21.305 

21.701 


Pressure ratio r = 1 + 


CIRCUIT LOSS FACTOR 



0.00 

0.004 

0.008 

0.0120 

0.0140 

0.0160 

0.0180 

0.0190 

0.020 

0.109 

0.121 

0.133 

0.153 

0.169 

0.186 

0.19f 


M = Test section Mach number 


Table 3. Tunnel and simulator test conditions. 




Fig. 

No. 

Total Temperature 

m 

Test Section 
(Mach Number) 

Total Pressure 
(atm) 

38 

100 

0.3 

1.57 

39 

100 

0.6 

1.57 

40 

100 

0.9 

1.57 

41 

100 

0.3 

3.00 

42 

100 

0.6 

3.00 

43 

100 

0.95 

3.07 

44 

100 

0.3 

5.00 

45 

100 

0.6 

5.00 

46 

100 

0.9 

4.93 

47 

200 

0.3 

1.57 

48 

200 

0.6 

1.50 

49 

200 

0.9 

1.57 

50 

200 

0.3 

2.65 

51 

200 

0.6 

3.00 

52 

200 

0.9 

3.00 

53 

200 

0.3 

5.00 

54 

200 

0.6 

5.00 

55 

200 

0.9 

5.00 

56 

275 

0 . 3 

1.57 

57 

275 

0.6 

1.57 

58 

275 

0.75 

1.57 

59 

275 

0.3 

3.00 

60 

274 

0.6 

3.00 

61 

275 

0.75 

3.00 

62 

275 

0.3 

5.00 

63 

275 

0.6 

5.00 

64 

275 

0.8 

5.00 
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Figure 3. Cooling capacity of liquid nitrogen. 
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p 

T 




m, 


'I P Tunnel average pressure, atm 

„ [ T Tunnel average temperature, K 

J — 1 1 


Plenum 




m. 


Fan 


Wg Mass Of tunnel gas. Kg 


Assumptions 

Perfect gas behavior 

Uniform tunnel temperature 

Work, potential & kinetic energy ignored 


first law of thermodynamics ... , 




+ u(m, - m„) 

L g 


Where 

h|^ Liquid nitrogen enthalpy J/kg 

u Specific internal energy 
hg Gas enthalpy J/kg 

Heat flow to tunnel metal wall J/sec 


Qp Heat flow from Fan operation J/sec 
ifig Gaseous nitrogen bleed kg/sec 

Liquid nitrogen mass flow kg/sec 


Figure 6. Thermodynamic model of a closed circuit 
cryogenic tunnel . 
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ENERGY STATE DIAGRAM 
CONTROL VEaOR MANIPULATION 



Mass — 


Figure 8. Tunnel state manipulation, temperature rise • 
at constant pressure. 


ENERGY STATE OiAGRAM 
CONTROL VEaOR MANIPULATION 

For a cryogenic pressure tunnel of fixed volume ideaiiy insulated from the ambient 

Rate control vectors 



Mass 


Figui'e 9. Tunnel state manipulation, temperature drop 
at constant pressure. 


ENERGY STATE DIAGRAM 
CONTROL VECTOR MANIPULATION 


For a cryogenic pressure tunnel of fixed volume ideally Insulated from the ambient 


Energy 



Mass 


Figure 10. Tunnel state manipul?.tion, pressure rise at 
at constant temperature. 
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ENERGY STATE DIAGRAM 
CONTROL VECTOR MANIPULATION 

For a cryogenic pressure tunnel of fixed volume Ideally Insulated from the ambient 

Rate control vectors 



Mass 


Figure 11. Tunnel state manipulation, pi’essure drop at 
constant temperature. 



T' Temperature K 


Figure 13. Enthalpy of aluminum alloy 6061 -T6 










Figure 15. 0.3-m cryogenic tunnel energy state diagram, 

ideal external insulation. 





Figure 16. NTF energy state diagram, ideal internal 
insulation. 
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Figure 18. Metal -to-gas heat transfer mechanism. 
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Figure 19. Spatial profile of total temperature and 
pressure. 


87 



CfT-*— 


'•mrv 




MACH NUMBER 
1.0 r 


MACH NUMBER PROFILE 
AREA RATIO BASED COMPUTATIONS 



M = i.O 




-NACELLE 


3rd 
CORNER 7 


4th 


TEST SECTION 1st 


#•1 


GN^ BLEED 


RNER 7 ^CORNER and PLENUM CORNER- r( 

SCREEN-' ^CONTRACTION '-IN JMIPmOf 


2nd 
■CORNER 


-CONTRACTION -LN^INJECTION^^^^ 




FAN 


DISTANCE ALONG CENTER LINE, m 


21.8m 


0.3 m TCT development 


Figure 20. Spatial profile of tunnel flow Mach number. 
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Figure 23. Spatial profile of local metal time constant. 
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Figure ^,4. Average metal time constant. 
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Figure 2b. Lumped multivariable model of a cryogenic 
tunne 1 . 
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Figure 27. Steady-state tunnel temperature gain for 
LNj input. 




Figure 29, 


Steady- state tunnel temperature gain for gas 
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Figure 30. 0.3-m --ransonic cryogenic tunnel control 

configuration. 
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Figure 51. Liquid nitrogen charge and gaseous bleed 
system schematic diagram. 






Operator display and control panel for hybrid 
simulation. 
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Figure 33. Hybrid simulation of a cryogenic tunnel. 
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Figure 34. Analog scheme for temperature d>Tiamics 
simulation. 
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Figure 36. Simulat 


or input function control 
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Figure 39. Tunnel and simulator responses; 0,6 M, 1.57 
atm, 100 K. 
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Figure 40. Tunnel and simulator responses; 0.9 M, 1.57 
atm, 100 K. 
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Figure 42 


Tunnel and simulator responses; 0.6 M, 5 atm 
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Figure 44. Tunnel and simulator respon 
atm, 100 K. 
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Figure 49. Tunnel and simulator responses; 0.9 M, 1.57 
atm, 200 K. 
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Figure 55. Tunnel and simulator responses; 0.9 M, 5 atm 
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Tunnel and simulator responses; 0.3 
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Figure 66. Warmup profile tunnel and simulator. 




